New experimental evidence that the proton develops asymptotically into a black disk 
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Recently, the Auger group has extracted the proton-air cross section from observations of air 
showers produced by cosmic ray protons (and nuclei) interacting in the atmosphere and converted 
it into measurements of the total and inelastic pp cross sections a to t and <7m e i at the super-LHC 
energy of 57 TeV. Their results reinforce our earlier conclusions that the proton becomes a black 
disk at asymptotic energies, a prediction reached on the basis of sub-LHC pp and pp measurements 
of a tot and p, the ratio of the real to the imaginary part of the forward scattering amplitude [M. M. 
Block and F. Halzen, Phys. Rev. Lett. 107, 212002 (2011)]. The same black disk description of the 
proton anticipated the values of cr tot and <Ti nc i measured by the TOTEM experiment at the LHC 
cms (center of mass) energy of yfs = 7 TeV, as well as those of o"i ne i measured by ALICE, ATLAS 
and CMS, as well as the ALICE measurement at 2.76 TeV. All data are consistent with a proton 
that is asymptotically a black disk of gluons: (i) both a to t and a- lnc i behave as In 2 s, saturating 
the Froissart bound, (ii) the forward scattering amplitude becomes pure imaginary (iii) the ratio 
finel/ftot = 0.509 ± 0.021, compatible with the black disk value of |, and (iv) proton interactions 
become flavor blind. 

PACS numbers: 12.38.Qk, 13.85Hd, 13.85Lg, 13.85Tp 
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Introduction. Recently, eight measurements of the pp 
total cross section have been made at energies beyond the 
Tevatron. At the LHC cms energy of 2.76 GeV, ALICE 
[l[ has measured <Ti ne u at the LHC cms energy of 7 TeV, 
ALICE [1], ATLAS 2], CMS @] and TOTEM Q have 
measured Cmci; crtot has been measured by TOTEM 
Most recently, the Pierre Auger Observatory has pub- 
lished pp cross sections for tr incl and a to t at 57 TeV 
The goal of this note is to update the evidence for the 
proton asymptotically becoming a black disk of gluons, 
using these new experimental results. 

At 2.76 TeV, the measured value of the pp cross section 

is: 

ALICE: cr incl = 62.1 ± 1.6(MonteCarlo) ± 4.3(lum.) mb. 

At 7 TeV, the measured values are: 

ALICE: dinci = 72.7 ± 1.1 (MonteCarlo) ± 5.1(lum.) mb, 

ATLAS: cr incl = 69.1±2.4(experim.)±6.9(extrapol.) mb, 

CMS: CTinoi = 68.0 ± 2.0(syst.) ± 2.4(lum.) ± 4(extrapol.) 
mb, 

TOTEM: cr ino i = 73.5 ± 0.6(stat.) ±* j(syst.) mb, 
TOTEM: cr tot = 98.3 ± 0.2(stat.) ± 2.8(syst.) mb. 

At 57 TeV, the measured values are: 
Auger: cr incl = 92 ± 7(stat.) ± 9 n (syst.) ± 7(Glauber) mb, 

Auger: cr tot = 133 ± 13(stat.)±2o(syst.) ± 16(Glauber) 
mb. 

We will show that all high energy measurements are 
in excellent agreement with recent predictions made by 



BH (Block and Halzen) |6j using a combination of an 
analyticity-constrained j7|j fit to sub-LHC data[l, Q to 
total cross sections and p- values in the energy range 6 < 
y/s < 1800 GeV, together with an eikonal model jlOl]. 

We will first recall the methodology used by BH for ob- 
taining accurate ultra-high energy extrapolations of the 
total and inelastic pp cross sections. 

Fitting data anchored by analyticity constraints — 
Using analyticity constraints [7] to anchor an analytic 
amplitude description of pp (pp) forward scattering am- 
plitudes [lH, BH Q made accurate predictions of the 
high energy behavior of both their total cross sections 
CTtot and p- values, 
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where /(#l) is the pp laboratory scattering amplitude, 
Ql, the laboratory scattering angle and p the labora- 
tory momentum. Saturation of the Froissart bound here 
means that the total cross section atot asymptotically 
behaves as In 2 s. Furthermore, the use of analyticity 
constraints allows one to anchor the fits to the accurate 
low energy cross section measurements between 4 and 6 
GeV, in the spirit of FESR (Finite Energy Sum Rules) 
7]. Making a local fit to the many very accurate exper- 
imental values of between 4 and 6 GeV, for both pp 
and pp, one obtains Q fixed "anchor-points" for a ± and 
their energy derivatives in Eq. (O at 6 GeV, the lowest 
energy value of our analytic amplitude fit. The model pa- 
rameterizes the even and odd (under crossing) total cross 
sections and p-values and fits 4 experimental quantities, 
o~p p (v), o~ pp (v), Pppiy) and p V p(y) to the high energy ana- 
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where the upper sign is for pp and the lower sign is for 
pp, and, for high energies, vjm ~ s/2m 2 . The even 
amplitude cross section is given by 



with v the laboratory energy of the incoming proton 
(anti-proton), to the proton mass, and the 'Regge inter- 
cept' p = 0.5. Note that vjm = s/2m 2 , so that asymp- 
totically, er° — ¥ In 2 s. 

The results Q of a global fit to both pp and pp p- values 
and total cross sections in the energy range 6 < ^fs < 
1800 GeV is shown in Fig. [TJ As seen from Eq. (j4j, 
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FIG. 2: Froissart-bounded analytic amplitude fits to the 
total cross section, a to t, for pp (dashed curve) and pp (dot- 
dashed curve) from Eq. Q, in mb vs. <*/s, the cms energy 
in GeV, taken from BH [SJ. The pp data used in the fit are 
the (red) circles and the pp data are the (blue) squares. The 
fitted data were anchored by values of of£ t and a^ t , together 
with the energy derivatives da^ t /dv and da^/dv at 6 GeV 
using FESR, as described in Ref. [9}. It should be noted 
that our ultra-high energy total cross section predictions that 
are made from our analytic amplitude fit use only total cross 
section data that are in the lower energy range 6 < \fs < 1800 
GeV. 
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FIG. 1: Froissart-bounded analytic amplitude fits to p, the 
ratio of the real to the imaginary portion of the forward scat- 
tering amplitude, vs. ^fs, the cms energy in GeV, taken from 
BH [9(. The pp data used in the fit are the (red) circles and 
the pp data are the (blue) squares. 

p — > as s — > 00, which is a requirement for a black disk 
at infinity. However, the tiny change in p from 0.135 at 
1800 GeV to 0.132 at 14000 GeV implies that we are 
nowhere near asymptopia, where p = 0. 

The fits for the pp and pp total cross sections are shown 
in Fig. [2j The dominant In 2 s term in the total cross sec- 
tion <t° (see Eq. ([S)) saturates the Froissart bound [l2j ]: 
thus it controls the asymptotic behavior of the cross sec- 
tions. 

Two low energy constraints on a^ t and er^, together 
with their energy derivatives da^/dv and da^/dv from 



FESR 9|) can be fixed precisely at 6 GeV by using the 
many accurate low energy total cross section measure- 
ments between of 4 and 6 GeV. These FESR con- 
straints can be used to fix the values of Co and (3-p> , two 
of the four parameters needed to determine a , the even 
high energy total cross section of Eq. (J5J) - These values 
of Co and f3-p> , together with with the 2 globally fitted val- 
ues of ci and ci required for <J°(v) (obtained from fitting 
the high energy total energy cross section and p mea- 
surements in the energy region 6 < y/s < 1800 GeV), are 
listed in Table HI We remind the reader that only data in 
the energy region 6 < ^/s < 1800 GeV are used in this 
global fit, together with the 4 to 6 GeV total cross section 
data used for the 6 GeV low energy 'anchor points'. We 
note that C2, the coefficient of ln 2 (s), is well-determined, 
having a statistical accuracy of ~ 2%. Thus, we see from 
Fig. [2J that the experimental data show that a saturated 
Froissart bound model is accurately satisfied for total 
cross sections Otot for both pp and for pp in the energy 
interval 6 < y/s < 1800 GeV; this accuracy of prediction 
mainly results from the use of the FESR constraints on 
the high energy analytic amplitude fit 0- 

Ultra-high energy cross total sections, for which there 
are no distinction between pp and pp interactions — both 
being given by <7° — are now completely predicted. For 
example, we obtain values for the total pp cross section 
of ct° = 95.4 ± 1.1 mb at 7 TeV [13] and 134.8 ± 1.5 mb 
at 57 TeV 14], compared to the experimental values of 



98.3 ± 2.9 mb and 133 ± 24 mb, respectively, where the 
systematic and statistical experimental errors have been 
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taken in quadrature. 



TABLE I: Values, in mb, of the parameters needed to deter- 
mine the even amplitude total cross section, <T°(y) of Eq. J5|, 
taken from Ref. [Sj]; for a fuller description, see the text. 



co=37.32 mb, 


/V =37.10 mb 


ci=-1.440± 0.070 mb, 


c 2 =0.2817± 0.0064 mb, 




Determination of the inelastic cross section — The in- 
elastic cross section, of nol , is determined by numerically 
multiplying the ratio of the inelastic to total cross sec- 
tion with the fitted total cross section a . The ratio of 
inelastic to total cross section was determined from an 
eikonal model, called the 'Aspen' model; for details see 
Ref. [10] and Ref. @|. 

This procedure is purely numerical; when fit by an 
analytic expression for the even amplitude high energy 
inelastic cross section of ncl (i/) given by 



I \ — flincl ( V \ , incl | inol , ( V \ 

^inelM = PV> ( — J + c + c l m ( — 

+4 ncl ln 2 (-), (6) 
\mJ 



we found that 
4^(u) = 62.59 (£) 
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24.09 + 0.1604 In 



+0.1433 In 2 (—) mb, 



(7) 



valid in the energy domain, y/s > 100 GeV. 

Results — The lower (red) plot of Fig. [3] is our pre- 
diction for high energy inelastic cross sections <7i, lc i as a 
function of the cms energy in the energy region 100 < 
y/s < 100000 GeV. The error bands corresponding to 
±1(7 are the lower dashed curves, and all of the existing 
inelastic cross section measurements for pp, as well as the 
5 new ultra-high pp measurements, are shown. Clearly, 
the agreement with experiment is excellent over the en- 
tire energy scale. 

Also shown as the upper (black) curve in Fig. [3] is 
our prediction for ultra-high energy total cross sections 
f tot (given by a of Eq. ([5]) using the coefficients of Table 
H|. Again, the error bands are the upper dashed curves. 
The excellent agreement with the new highest energy (57 
TeV) experimental measurements is striking. Since none 
of the experimental datum points in Fig. [3] are used in 
making these predictions, it is clear that our In 2 s pre- 
dictions for (Tinei and a tot are strongly supported by the 
existing ultra-high energy measurements. Further confir- 
mation (at 14 TeV) is possible in several years, when the 
LHC runs at its design energy. 



FIG. 3: Predictions for crtot and <7i n ei vs. y/s for pp and 
pp. For crtot, we have compared our predictions with recent 
pp TOTEM data at 7 TeV and Auger data at 57 TeV, while 
for (Jinoi, we have compared our results with 2.76 TeV pp data 
from ALICE, 7 TeV pp data from ALICE, ATLAS, CMS and 
Totem, as well as with the 57 TeV pp inelastic cross section. 
The upper solid (black) curve is the central-value prediction 
for a t t and the lower solid (red) curve is the central-value 
prediction for (Tinci- The dotted curves are the errors (±1<t) 
in our predictions, due to the correlated errors of the fitting 
parameters. We emphasize that none of the datum points in 
this plot have been used in our predictions. 



Asymptopia — Finally, we Q determined the ratio of 
^incl (s)/otot(s) as s 
coefficients in crP nol and a 



oo , given by the ratio of the In s 
respectively, i.e., 



Omcl 
ftot 



1433 

^— = - = 0.509 + 0.021, ass->oo,(8) 

c 2 0.2817 ' ' w 



that is well within error of the expected value of \ that 
is appropriate for a black disk at infinity. 

These new experiments confirm our earlier results Q 
that the asymptotic proton is made up of gluons and thus 
is flavor blind, yielding the same asymptotic cross sec- 
tions for pp, Kp, pp, irp and, through vector dominance, 
to 7p and j*p scattering, as in deep inelastic scattering. 
We also found [(| that the coefficient C2 of the In 2 s term 
in the total cross section corresponded to a lowest lying 
glueball mass of M glueb all = (2tt/c 2 ) 1/2 = 2.97 ± 0.03 
GeV. 

The experimental ratio of cxinei/ctot ~ 0.72 at 7 TeV; 
at the highest available energy of 57 TeV, the ratio very 
slowly decreases to « 0.69, not even close to the asymp- 
totic limit of 0.5. Hence, even at the highest cosmic ray 
energies, we still are a very long way from asymptopia 
and will never get much closer to it. However, it is most 
interesting that the essential principles of analyticity and 
unitarity, the underpinnings of our theoretical results, arc 
experimentally validated up to these ultra-high energies. 
Clearly, it will become the task of lattice QCD to ex- 
tend these results to the enormous energies needed to 
approach asymptopia. 
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Additional experimental confirmation that the pro- 
ton asymptotically approaches a black disk is given by 
Shegelsky and Ryskin [15J , who analyze all available high 
energy data on the energy dependence of the shrink- 
ing of the diffraction cone (the nuclear slope parameter 
B = <icr c i as tic In £ | t = o) i finding agreement with our re- 
sults. 

Conclusions — We find that: 

1. both the total pp cross section and the inelastic 
cross section are fit up to y/s — 57 TeV by a sat- 
urated Froissart-bounded In 2 s behavior that is as- 
sociated with a black disk. 

2. the forward scattering amplitude is pure imaginary 
as s — ¥ oo, as is required for a black disk. 

3. the ratio of o'mci/o'tot — » 0.509±0.021, as s — > oo, 
compatible with the black disk value of 0.5. 

Thus, we conclude that existing experimental evidence 
strongly supports the conclusion that the proton becomes 
a black disk at infinity. Our result may have implications 
for theories of "new" physics; for instance, string theories 
whose additional dimensions modify forward scattering 
amplitudes as shown originally by Amati, Ciafaloni and 
Veneziano 16]. 
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